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. Abstract 

\ In this article we considered the integrable problems of three vortices on 

a plane and sphere for noncompact case. We investigated explicitly the prob- 



<^ . lems of a collapse and a scattering of vortices and obtained the conditions 



of its realization. We completed the bifurcation analysis and investigated for 
collinear and Thomson's configurations the dependence of stability in linear 
approximation and frequency of rotation in relative coordinates from value of 
pH \ a full moment. We indicated the geometric interpretation for characteristic 

situations. We constructed a phase portrait and geometric projection for the 
integrable configuration of four vortices on a plane. 



1 Motion of vortices on a plane 

Let us consider motion of vortices on a plane under the following condition [at are 
inverse intensities of point vortices = l/Ffc) 

A = Aia2 + 02^3 + "lOs < . (1) 

In case A < algebra of Poisson brackets is the algebra M © so(2, 1), and in case 
A = it is solvable algebra. As indicated in ^2]; the symplectic sheet in both cases 
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is noncompact (in the first case the symplectic sheet is hyperboloid, in second it is 
paraboloid), and the trajectories of a representing point on it can be both finite and 
infinite (this also concerns the dynamics of three vortices in relative distances). In 
the second case one can speak about a scattering of vortices. 

Let us consider, at first, the case ^4 < 0. Without loss of generality let us put 
Fa > 0, Fi < 0, -Fi < F2 + F3 > 0. Further, D is the central hnear function 
of algebras of Lie-Poisson brackets (D — Y^a^M]^, are quadrates of mutual 
distances between vortices). From 

2 —8010203 3 
A = (01 + 02 + 03), (2) 

for a stability coefficient of Thomson's configurations, it follows, that the stabihty of 
such configuration in a linear approximation for various D is determined by sign of 
the sum of inverse intensities S = XI '^fc- Therefore, we shall consider three separate 
cases corresponding to the values S > 0, S < 0, and S = 0. 

The geometric interpretation and the bifurcation diagrams for various charac- 
teristic combinations of parameters A, S and D are indicated correspondingly in a 
Fig. 1 and Fig. 2. 

Case A < and S > 0. The geometric interpretation shows that in this case for 
D < only finite motions are possible for which trajectory on a plane of an integral 
of the full moment is limited from both sides by the condition > 0. 

The bifurcation diagram (Fig. 2a) in this case reminds the diagram for vortices 
on a plane in compact case for positive intensities with the exclusion, that there is 
only one collinear configuration. Another difference from compact case is that the 
type of stability varies: Thomson's configurations are stable for positive intensities in 
compact case and are unstable in this case (Fig. 3a), and the collinear configurations 
are stable. Both Thomson's, and the coUinear solutions are determined only in area 
D > (Fig. Ic). For positive values of the full moment both finite, and infinite 
motions are possible. By virtue of instability of Thomson's solution, the solutions 
with any positive energy are allowed. 

Case A <0 and S < 0. In this case for any value of the full moment the motion 
is finite (Fig. Id-f). Thomson's and the coUinear solutions are determined only for 
D < 0. That solutions are remarkable because the energy of all configurations tends 
to infinity with decrease of the absolute value of the full moment down to zero (see 
Fig. 2b). The curve of Thomson's configurations is situated above of the curve for 
collinear configurations and limits from above the area of possible movements. They 
are accordingly stable and unstable in a linear approximation (Fig. 3b). 

Case A < and S — 0. For the value D < the tangency of the trajectory 
with the boundary of area is possible. This trajectory corresponds to stable (see 
Fig. 3c) collinear configuration (Fig. 2c), which energy tends to zero with decrease 
of D. Thomson's solutions are presented only ior D — and are degenerated, and 
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have a neutral equilibrium (A^ = 0). Such solutions are possible for any distances 
between vortices and fill the whole straight line at Fig. Ih. The dependence of an 
angular velocity from the distance between vortices is determined by the formula 



where M is quadrate of the distance between vortices. This straight line, corre- 
sponding to Thomson's configurations, divides areas of scattering and collapsing 
behavior of three vortices. Under a "collapse" we understand the process of simul- 
taneous collision of (in this case, three) vortices. More detailed study of emergence 
of a collapse is indicated further. 

For D 7^ the finite motions are presented for any value of the moment, and 
infinite motions occurred only for negative values. 

The dependence of an angular velocity of rotation both Thomson's and coUinear 
configurations in first and second cases is qualitatively identical and is defined by 
functions monotonically decreasing, with increase of absolute value of a moment. 
For case S* = the angular velocity of rotation of a collinear configuration is slowly 
increased with increase of the absolute value D. 

The last case corresponding to noncompact motion of vortices arises under con- 
dition of A = (solvable algebra). The motion is possible only for positive values of 
D. There are only two stationary configurations: one Thomson's and one collinear 
(Fig. 2d). For D = Q three vortices are placed on one straight line, and each vor- 
tex is placed at center of a vorticity of two remaining. It is a corollary of general 
possibility of reducing of a problem of n -f 1 vortices to problem of n vortices. 

2 Motion of vortices on a sphere 

Similarly to Section 1 we shall construct the bifurcation diagrams for the sphere, and 
also plots of absolute motion, using for small D appropriate dependencies for a plane 
and the method of a continuation on the parameter. The values dk = 4:R^{ai + aj) 
determine those values of the full moment, which can arise in problem of two vortices. 
The geometric interpretation for motion on the sphere is indicated in a Fig. 5. 

Case A < and S > 0. In addition to Thomson's and collinear configurations 
that exists in plane case, for D > there are two collinear configurations, appearing 
from a problem of two vortices (Fig. 4a). By continuation on the moment one of them 
merges with Thomson's, which disappears for value dr- The collinear configuration 
continues and merges with another collinear configuration, that has an analog on 
the plane. After this both configurations disappears also. For D > d^ and D < di 
the stationary configurations do not exist. Thomson's configuration are in this case 
always unstable (Fig. 6a), and the collinear configurations are stable. 




(3) 
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Case A < Q and S < 0. In this case the confluence of Thomson's and one of 
coflinear (closest to Thomson's on energy) configurations (Fig. 4b) arc happened. 
Besides there are colhnear configurations from a problem of two vortices (for one 
of them £^ I— > oo by coincidence of vortices because of existence of one negative 
intensity). One of them extended into area of positive values of D. Thomson's 
configuration, as well as in plane case, is stable in the linear approximation (Fig. 6b). 
Some of collinear configurations are also stable. One of them becomes stable in 
confiuence with Thomson's configuration, and other loses stability in confluence 
with unstable collinear configuration. 

Case A < and S — 0. Instead of one collinear configuration in case of a plane 
in this case there are three various branches of colhnear configurations, as for D > 
and for D < (see Fig. 4c). The solution with higher energy is unstable for D > 
and disappears merging with collinear configuration with lower energy. The solution 
with the lowest energy is stable. For D < situation is reversed. Thus collinear 
solution with lower energy merges with the solution, appeared at D — d2- As well 
as on the plane, the Tomson solution are present only ior D — 0. In this case all 
trajectories are collapsing (Fig. 5h). 

At last, under condition of A = 0, the bifurcation diagram is indicated in a 
Fig. 4d. As well as in plane case, the stationary solutions exist only for D > 0, 
(Fig. 5 j-1). The difference consists in emerging of collinear configurations of a 
problem of two vortices. One such conflguration merges with Thomson's, then both 
disappear. All collinear solutions are stable, and Thomson's is unstable (Fig. 6d). 

The dependence of angular velocity of rotation for the listed situations is repre- 
sented in a Fig. 7. As well as in compact case, the frequency of rotation increasing 
with creation of the third vortex from the problem of two vortices. 

Remark 1. Such magnification, as well as in compact case, is explained by the 
fact, that for the birth of a pair of vortices from one vortex of total intensity there is 
a rotation of a collinear configuration around an axis, which is passing through the 
third vortex. In the result this vortex actually does not influence on rotation of two 
appearing vortices, which angular velocity depends on the mutual distance between 
them as 



Where M3 is quadrate of the distance between first two vortices, and ri,r2 are 
intensities, corresponding to it. As the last expression shows the angular velocity 
tends to infinity with M 1— > 0. 

The study of dynamics of angles of dechnation of rotation axis to the plane of 
vortices shows, that in first and second listed cases the angle of declination of axes 
increasing from zero up to 71 /2. So, therefore the rotation axis tends to take a 
position in the plane of vortices and appropriate Thomson's a configuration turned 
into collinear (Fig. 8). 
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Remark 2. Note, that a condition of existence of static configurations on a 
sphere for a specific collection of intensities is 



010203(01 + 02 + 03) > , 



(5) 



which determines stability of Thomson's configuration 



(oi + 02 + O3) 



010203(01 + 02 + 03) . 



(6) 
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This correlation is easy for understanding from the geometric interpretation (see 
Fig. 5d), which illustrates a possibility of reaching of function H to its extremum. 
Let us note, that as well as in compact case, that was considered in 11 g) /. the static 
configurations on a sphere are observed only under fulfillment of the condition rir2 + 



r2r3 + riFg > o, (Ftg. 7). 

3 The condition of a collapse of vortices on a 
plane and a sphere 



Let us consider conditions of emergence of a collapse of three vortices on a plane and 
a sphere. It is known, that the collapse is impossible for any pairs of vortices taken 
separately, because with the approach of such pair influence of the third (remote) 
vortex negligible small, and two vortices on a plane (sphere) move relatively each 
other so, that the distance between them did not change. The necessary condition of 
the collapse is the fulfillment of D = 0, since it follows Mk = 0. The condition D = 
allows to proceed from four- dimensional to three-dimensional Lie algebra, for which 
actual dynamics takes place on the singular symplectic sheet, so that the nontrivial 
singular symplectic sheet (surface of a cone) will correspond only to algebra so(2, 1). 
The similar statements are valid also for a simultaneous collapse n vortices, which 
is not almost investigated (except for cases = 4, 5) [TJ[TU]. 

The conditions of a collapse of vortices on a sphere, obviously, will coincide 
with conditions for a plane, since on small distances the influence of curvature to 
dynamics of vortices is insignificant. 

Remark 3. The problem of the collapse is one of the most interesting problems 
connected to vortices and represents large interest for the theoretical hydromechanics 
as one of the models, which can be used for understanding of the transition to the 
turbulence contained in a nonuniqueness of solutions of the hydrodynamic equations 
of the Euler. Really, the theorem of existence and uniquenesses for these equations 
are proved in the supposition of sufficient smoothness of an initial field of velocities. 
From a mathematical point of view the process of the collapse of vortices represent- 
ing confluence of the special solutions of Euler equation of 5 -function type, under 
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inverting of time will determine disintegration of vortices with corresponding loss of 
uniqueness. Therefore large interest represents the study of this problem from the 
point of view of a regularization of collisions in the similar way, how it is done in 
classical celestial mechanics For physics of atmosphere the phenomenon of the 
collapse can be considered as a model of forming of large atmospheric vortices. 

We research at first possibility of a liomogeneous collapse for three vortices. 
For the homogeneous collapse all distances between vortices depend on time in the 
identical way and are in the constant proportion ITUj. Using isomorphism with a 
problem of Lottka-Voltera (both for a plane, and a sphere jHI), we shall consider a 
homogeneous system of equations in the form 

Ml = riMi(M2 - M3) . (7) 

The homogeneous asymptotic solutions ((Tj) we shall search as 

Mk = —. (8) 
r 

Let's note, that such solutions are used in Kovalevskaya method for construction 
of full-parametric Laurent expansion. The nontrivial solutions of such form are 
possible only under condition of 5* = 0. If this condition is fuUfiled, the solution can 
be written as 

C ,^C + a3 ^ ^ C — 0,2 /X 
Mi = -, M2 = M3 = 9 

T T T 

where C is arbirary constant. 

The solution (jH)) is full-parametric and contains only collapsing and running 
up (for a plane) trajectories. However, there is still the special solution Mi = 
M2 = M3 = const, describing the Thomson's configuration, which in this case are 
degenerate. Using the relation 

dt = ^ ^ dr, (10) 

it is possible to receive an asymptotics of the solution Q in real time 

1. for a plane: t = ^, Mk = C^t; 

2. for a sphere: t = AR^ (l - ^1-^), M, = D^R^ (l - (l - ^)'), 

where A, B, Cl, D, = const, R is radius of a sphere. For a sphere the absolute 
motion of vortices for conditions D = and S = consists in scattering of vortices 
from one point up to the moment of reaching of equator and further approach in 
some other point. Let us analyze the possibility of the collapse in the system of three 
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vortices on a plane in general inhomogeneous case. Writing the necessary condition 
of the collapse D = in absolute variables 

'J2rk)i-P'-Q' = 0, (11) 

where the uninvolute integrals of motion P,Q,I are equal 

P = J2^kXk, Q = J2^kyk, I = J2^kixl + yl), (12) 

we found two special cases ^ = and 7^ 0. Without loss of generality let 

us assume Fi = —1, F2,F3 > 0. Under conditions ^Fj^ = 0, from (fT^ it follows 
that P = Q = 0. It means, that the third vortex is at center of vorticity of first two 
vortices, rotating uniformly with frequency 

1 Fi + F2 / Fi F2 
27r Ms V Fi 

around the point with the radius- vector 
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where M3 is quadrate of the distance between the first and second vortices, and 
fk, Ffc are the radius-vector and intensity correspondingly. Therefore, collapse for 
case of zero total intensity is impossible. 

Under condition of ^ F^ 7^ the collapse is possible only under condition of 
noncompactness of algebra of vortices A = 0102 — ai — 02 < 0, as in compact case 
a symplectic sheet (the sphere) does not pass through a beginning of coordinates 
Mk = 0. 

For the determination of sufficient conditions of the collapse we shall consider 
the projection of trajectories on a plane Mi, M2. Let us express M3 from the integral 
of the full moment (for D = 0) 

Ms = aiMi + a2M2 . (13) 

The physical area on a plane Mi, M2 for D = is limited by straight lines, passing 
through the zero with coefficient of declination 

A,,, . (IMEI) . ,14) 



1 - a2 

The trajectory of regularized system is determined by the relation (|T!^ and equation 
for energy of the system, which can be presented as 

Mf''iM2"''^(aiMi + 02 M2) = const (E) . (15) 

The analysis of the equation (fT3j) shows, that for various values of parameters ai, 02 
there are three various kinds of trajectories: 
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1. if Oi + 02 < 1, then all trajectories are compact, looks like closed loops, 
which are going out the origin of coordinates, tangentcing axes OMi and OM2 
(Fig. 9a, curves b, c); 

2. if Oi + 02 > 1, then all trajectories are noncompact, going to infinity, asymp- 
totically press oneself to axeses OMi and OM2; 

3. if Oi + 02 = 1, then all trajectories are straight lines which are going out from 
the origin of coordinates under different angles (the Fig. 9a, curve a); 

Let us mark on a plane of parameters Oi, 02, areas appropriate to given types of 
trajectories. It is necessary to mark on this plane also values of parameters describing 
various forms of physical area. By (fT^ . the physical area for oi 7^ 1, 02 7^ 1 (such, 
that A < 0), represents an interior of an sharp angle located inside a quadrant 
Ml > 0,M2 > 0. If ai = 1 (02 = 1), one of the sides of physical area coincides 
with an axes OM2 {OMi). In case ai = 02 = 1 the motion is allowed in the whole 
quadrant Mi > 0, M2 > 0. 

Comparing possible types of trajectories (for A < 0, D = 0) with types of areas 
of motion with consideration, that on reaching of the boundary motion continues 
on the same trajectory in the opposite direction, we conclude: 

1. In the system of three vortices the homogeneous collapse (scattering) is possi- 
ble, it happens when the relation ai + 02 = 1 is executed for inverse intensities; 

2. Scattering of vortices is possible only under conditions ai = 1, 02 = 1, ai = 
02 = 1 (in this case the vortex motion never passes through the coUinear 
configuration) ; 

3. For other values of inverse intensities the motion of vortices is bounded 
between two collinear configurations, and the distance between them is limited 
(see Fig. 9b). 



4 Scattering of vortices on a plane 

In the system of three vortices we shall name scattering trajectories, such trajec- 
tories that at least one of mutual distances (Mi,M2,M3) is infinitely increased. In 
difference from the collapse, the scattering can happen for D ^ Q. 

Let us define new variables, with the help of which scattering problem of vortices 
on a plane can be reduced to the research of collapse: 
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Their equations of motion have a form 

^ xi = xi{{T2 - T3)xi + T3X2 - T2X3) , (17) 

and the bounding physical area relation > 0, will not be changed 

2{xiX2 + X1X3 + X2X3) — xl — xl — xl> . (18) 

Let us select, as above 03 = —1, ai,a2 > 0. The trajectory in space of variables 
Xi,X2,X3 is set by integrals of the moment D and the energy E, which can be 
presented in the form 

aiXi + a2X2 + a^xs = D^x\X2Xz , (19) 



^a2-ai-l^ai-a2-l^l+ai+a2 _ ^^^^ ^ (^20) 

where CiE) is some function of energy. 

For D = the type of trajectories, defined by (fTIH) . is similar to the one investi- 
gated above. In this case to scattering 2) in a system of three vortices corresponds an 
inhomogeneous collapse in a system (fTTjl . for conditions ai = 1, 02 = 1, ai = 02 = 1 
(see Fig. 9b). 

It is interesting to note, that the homogeneous collapse and the scattering in 
variables remains homogeneous in variables x^, only direction of motion on 
trajectories varies. 

Obvious quadratures and the analysis of absolute motion is in this case indicated 
in 1^. 

The numerical research for 7^ shows, that in variables there are only 
trajectories of types 1), 3) of the previous section. In this connection, it seems, 
that the availability of vortex pairs {a\ = 1, 02 = 1, ai = 02 = 1) in the system is 
necessary and sufficient condition of scattering. 

Remark 4. In a work 16] reducing of the order is executed and the phase portraits 
of integrable case of a problem with zero total intensity and zero total moment are 
indicated. In this work the canonical form of equations of motion is used. For 
the indicated conditions on a level of algebra of brackets there is an reduction of a 
problem of four vortices to a problem of three vortices with the reduced Hamiltonian. 
On a symplectic sheet of three-dimensional vortex algebra, defined by introduced 
(as a corollary of the integrability) an invariant relation such as an integral of the 
moment with a constant Di, standard canonical variables L,l are presented. The 
condition of compactness of a symplectic sheet will have a very simple form, if three 
vortices of four have intensities of same sign. In compact case, for various values 
of intensities the phase portraits are represented at Fig. 10. In a general situation 
of unequal intensities there are six collinear stable solutions and three uncollinear 
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unstable ones. The last solutions generalize of Thomson's configuration, however 
distances between vortices are not equal. The connection between energy E and 
moment Di for these solutions is determined by relation 

where /(Fi, Fa) is some function dependent on intensities. The bifurcation analy- 
sis, consisting in a determination of an explicit form of function /(Fi, F2, F3), can be 
executed similarly to Section 1. In case if there are only two intensities of the same 
sign, the symplectic sheet is noncompact and the scattering is possible. The regular 
scattering, for example, is possible in case of interaction two (generally speaking, 
various) vortex pairs As far as we know, in general case the conditions of a 
collapse and a scattering in considered problem are not investigated. 

After submission of the article in the journal the authors have received prints of 
articles [1^1 ITil [T3] from J.E.Marsden and P.K.Newton. Their results were received 
with the help of the "canonical" approach and are contained some results of our 
work. We are extremely grateful to them. The authors thank I. S. Mamaev and 
N. N. Simakov for useful discussions and help in work. The work is carried out 
under the support of Russian Fond of Fundamental Research (96-01-00747) and 
Federal program "States Support of Integration High Education and Fundamental 
Science" (project No. 294). 
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Figure 1: The geometric interpretation for various values of parameters A, S, D. 
The dark color designates area of positive values > 0, for which > 0. For 
figure h) the bold line corresponds to a set of Thomson's configurations. 
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Figure 2: Bifurcation curves on a plane for cases a) A < and > 0; b) A < and 
S* < 0; c) A < and S' = 0; d) A = and S > 0. Here and further dot-dashed hne 
corresponds to Thomson's configurations, and sohd hues to colhnear. 
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Figure 3: The stabihty coefficient on a plane for cases a) A < and 5" > 0; b) A < 
and 5 < 0; c) A < and 5 = 0; d) A = and S > 0. 
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Figure 4: Bifurcation curves on a sphere for cases a) A < and S > 0; h) A < 
and 5 < 0; c) ^ < and 5 = 0; b) ^ = and S > 0. 
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Figure 6: Stability coefficient on a spfiere for cases a) A < and S > 0; h) A < 
and < 0; c) A < and = 0; b) A = and S > 0. 




Figure 7: An angular velocity on a sphere for cases a.) A < and S > 0; h) A < 
and 5 < 0; c) A < and 5 = 0; b) A = and S > 0. 
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Figure 8: The angle of declination on a sphere for cases a) ^4 < and 5" > 0; b) 
^ < and 5 < 0. 




Figure 9: 1) Possible asymptoticses of a behavior of trajectories of the system of three 

vortices near to zero on a plane (Mi,M2) in case of the collapse (the nonphysical 
area is shaded); 2) Planes of parameters (ai, 02 for 03 = 1). The collapse is possible 
on straight line 02 = ai — 1, scattering corresponds to asymptotes ai = 1, a2 = 1. 




Figure 10: The geometric projection for case of four vortices on a plane with zero 
full intensity for a) oi = —02 = as = —04; b) oi = 02 = 03 7^ 04; c) ai 7^ 02 7^ 03. 
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Figure 11: The phase portrait of motion of four vortices on a plane with zero total 
intensity for a) al ^ a2 ^ a?) ^ al] b) al = a2 — a3. 
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